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Abstract 

We propose a classical analogue of the vertex algebra in the context of classical inte- 
grable field theories. We use this fundamental notion to describe the auxiliary function 
of the linear auxiliary problem as a classical vertex operator. Then using the underlying 
algebra satisfied by the auxiliary function together with the linear auxiliary problem we 
identify the local integrals of motion, which by construction are in involution. The time 
components of the Lax pair are also identified in terms of the classical vertex operators. 
Systems in the presence of point like defects as well as systems on the semi-infinite line are 
investigated. Specific examples associated to the classical Yangian and twisted Yangian 
are also presented. 
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1 Introduction 

The starting point in the present analysis is the notion of vertex operators at the quantum 
level and the associated Faddeev-Zamolodchikov algebra [1], Inspired by these ideas we 
propose a similar algebraic formulation to deal with classical integrable field theories on 
the infinite or semi-infinite line. It is worth noting that such ideas at the classical level 
were briefly discussed in [2], but a systematic construction of classical vertex operators, 
the generating function of the local integrals of motion as well as a construction of the 
time component of the Lax pair in terms of the classical vertex operators was not really 
demonstrated. We should stress that one of the key points of the present analysis is 
the identification of the auxiliary function of the auxiliary linear problem as the classical 
version of the vertex operator. Moreover, particular emphasis in the present investigation 
is given in integrable systems in the presence of point like defects as well as in the presence 
of non-trivial boundary conditions. 

Recall that vertex operators <f> at the quantum level satisfy the Faddeev-Zamolodchikov 
algebra [1] 

<5 < i2(Ai — A 2 ) <hi(Ai) $ 2 (A 2 ) = 4*2(A 2 ) <hi(Ai), (1-1) 

where S is the physical scattering matrix S £ End(C- / ^(8)C A ^), solution of the Yang-Baxter 
equation, and $ is the quantum vertex operator which is an W-dimensional vector with 
elements being realizations of the underlying affine quantum algebra. Here we are going 
to formulate a classical analogue of the latter quantum objects <f> as well as the algebra 


1 






(1.1). This will be basically achieved by means of certain vectors -tensor products of 
which are eigenvectors of the classical r-matrix- suitably implemented at the ends of the 
system. 

The paper is organized as follows: in the next section we briefly review the auxiliary 
linear problem and also recall the notion of the classical monodromy matrix together with 
the corresponding classical quadratic algebra. Then based on these ideas we introduce the 
classical analogue of the vertex operator, which coincides essentially with the auxiliary 
function of the auxiliary linear problem, which in turn satisfies a classical version of 
the vertex algebra (1.1). We then introduce the generating function of the integrals of 
motion in terms of classical vertex operators. The extracted charges are in involution by 
construction due to the fact that the classical vertex operators satisfy a quadratic vertex 
algebra. Based on the classical vertex algebra and the Hamiltonian equations of motion 
we derive the time components of the Lax pairs. This setting is generalized for systems 
in the presence if point like defects as well as for systems on the half line. In section 3 
we exemplify our construction using as a paradigm a whole class of systems associated to 
the classical Yangian. Explicit expressions of the generating functions of the integrals of 
motion as well as the time components of the Lax pairs are provided. 


2 The general setting 

The starting point in our analysis will be the auxiliary linear problem, consisting of the 
Lax pair U, V and the pair of differential equations (see [3] and references therein): 

— = U(x,t;A)4- 

— =Y(x,t,\) T. (2.1) 

U, V are in general AT x AT matrices with entries being functions of x, t and depending 
also on a spectral parameter. The main proposition of the present investigation is that 
the auxiliary vector functions T are essentially classical variants of the vertex operators, 
and they satisfy a classical version of the vertex (Faddeev-Zamolodvhikov) algebra (1.1). 

Let us recall the classical monodromy matrix, which will be an essential object in our 
formulation 

T(x, y, t; A) = V exp j J dx' U(x / , t\ A)j, (2.2) 

is a solution of the first of the equations (2.1). Moreover, the monodromy matrix satisfies 
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the quadratic classical algebra 1 

{Ti( A,), T 2 (A 2 )} 


ri 2 (A!-A 2 ), Ti(Xi) T 2 (A 2 ) , 


(2.3) 


where as usual we introduce the notation: Tj = T<g)I, T 2 = I(g)T and so on. The classical 
r-matrix satisfies the classical Yang-Baxter equation [4] 


r i 2 (Ai — A 2 ), r 13 (Ai) + r 23 (A 2 ) + r 13 (Ai), r 23 (A 2 ) 


= 0 . 


(2.4) 


We shall be mainly interested henceforth in integrable models on the full line or on 
the semi-infinite line. We shall begin our derivation with integrable systems on the full 
line considering the interval [-L, L\. We define the following basic objects 

T(x, —L, t] A) = T(x, —L, t: A) V 

T*(L, x, t; A) = V* T(L, x, t; A), (2.5) 

where V, V* are Af row and column vectors respectively such that: 

V* V 2 * r 12 (A) = /(A) V* V*, r 12 (A) W Y 2 = /(A) W Y 2 (2.6) 

/(A) a function of the spectral parameter depending on the choice of the r-matrix. 

Then it is straightforward to show using the algebraic relation (2.3) and the definitions 
(2.5), (2.6): 

— ti2(v - v) *i'i(41j 4 / 2 ( 42 ) 

{'U(V), 4 * 2 ( 42 ) } = -*;(4 1 ) 4 * 2 ( 42 ) 1,2(4! - 4 2 ) (2.7) 

(see also similar relations at the classical level in [2]), we define the “shifted” classical 
r-matrix: 

t(A) = r(A) - /(A) (2.8) 

Note that we consdider here a quite generic classical r-matrix: r G End(C Ar <g) C^). In 
the last section we are going to restrict our attention to the classical Yangian. 

Keeping in mind the latter relations we can introduce the fundamental quantity: 

T(L, —L, t; A) = 0, t A) T(0, -L, t A). (2.9) 

The latter (2.9) provides the charges in involution of the system under consideration, 
indeed it may shown via (2.7) that 

{Tx(A 2 ), T 2 (A 2 )} =0. (2.10) 

1 Note that throughout the text the dependence on x , t is implied even if it is not explicitly stated. 
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Expansion of the ln(TT) clearly provides the local charges in involution (local integrals of 
motion). Later in the text we are going to consider a class of integrable models associated 
to classical Yangians [5], with typical examples e.g the generalized NLS and the Landau- 
Lifshitz models. 


Employing the formulation described above we may extract the hierarchy of charges as 
in the inverse scattering process. The main advantage here is the fact that the auxiliary 
function may be thought of as a classical vertex operator satisfying the classical algebras 
(2.7), and thus involution of the charges is by construction guaranteed. Moreover, due to 
the fact that the classical vertex operators are immediately associated to the monodromy 
matrix facilitates technically the derivation of the desired quantities. The time depended 
equation of the pair (2.1) is not really used as opposed to the inverse scattering method¬ 
ology where the time dependent equation is exploited in order to prove the conservation 
of the extracted local charges. In fact, the time components of the Lax pair associated to 
each integral of motion can be also derived in terms the classical vertex operators. 

Indeed, we shall derive below the V-operators in terms of the classical vertex operators 
exploiting the underlying algebra. Recall the fundamental Poisson commutator (see e.g. 

[ 3 ]) 


|Ti(t4, B,t) A), U 2 (x,t; 
x e [A, B], 


dM 12 (x,t] A,/i) 

dx 


Mi 2 (x,t; A,/i), 


U 2 (sc,t;/i) , 

( 2 . 11 ) 


where we define 


Mi 2 (x,i; A, /j) = Ti(A,x,t\ A) r n {\- /() Ti{x, B,t~, fi) 


( 2 . 12 ) 


Recall also the definitions for T and Y* and formulate: 

<9N 2 (x, t ; A, //) 


n(L,0,t-X) ^(0,-L,f;A), U 2 (x,t; 
where we define 


dx 


+ 


N 2 (x, t; A, /i), U 2 (x,t;/i) 
(2.13) 


N 2 (x, t: A, p) = ^l(L, x, t; A) r i2 (A - /i) Ti(x, -L, U A) 


(2.14) 


We wish to consider ln(T) as the generating function of the local integrals of motion 
as well as the zero curvature condition arising as a compatibility condition of (2.1) 


V\x, t; A) — U(x, A) + Y(x,f;A), U(x,f;A) 


= 0 . 


(2.15) 
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the “prime” denotes derivation with respect to x, and the “dot” denotes derivation with 
respect to t. Then 


{in 0, t] A) Ti(0,-L,t; A)), U 2 (x, t; /i)| = 

where 


dV 2 (x,t, A,/z) 
dx 


V 2 (x,t; A,/i), HJ 2 (x,t;h) 
(2.16) 


V 2 {x,t- \,n) = T(A) 1 A) ri 2 (A - fi) Ti(a;, -L,t] A) (2.17) 


Expansion of V in powers of A 1 provides the hierarchy of V-operators associated to each 
integral of motion. 


2.1 Defects 


The approach described above may be easily generalized in the presence of local integrable 
defects (see also e.g. [6]—[9] for recent results in this context). In fact, the relevant object 
in this case, i.e. the generating function of the integrals of motion would be 

©(A) = T* + (L, x 0 , t ; A) L(x 0 , t; A) T _ (xo, —L, t; A) (2-18) 


the plus and minus superscripts in the expression above refer to the left right bulk theories, 
the defect matrix L is of the generic form 

a r 

L(A) = ^L^(A) ea, (2.19) 

*.j=i 


Cij are in general Af x Af matrices such that: (e^)*,* = So- 8ji . L satishes the quadratic 
algebra 


Li(A 2 ), L 2 (A 2 ) 


ti2(Ai — A 2 ), Li(A 2 ) L 2 (A 2 ) 


( 2 . 20 ) 


Exploiting the algebraic relations between the vertices (2.7) as well as the ones regarding 
the defect matrix (2.20) we conclude as expected that © provides the charges in involution, 


{b(A0, D(A 2 )} =0, (2.21) 

i.e. the involution of the charges derived via the expansion of the generating function D 
is guaranteed by construction. 

We shall derive below the V-operator associated to the defect point computed from 
the left and from the right in terms of the classical vertex operators. The V-operators at 
x ^ Xq coincide with the bulk expressions for the left and right theories, and we do not 
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rederive them here for brevity. To derive the V-operator on the defect point we recall the 
zero curvature condition at x = x 0 [8] 

d L( y A) = ¥ + (x 0 , t- A) L(x 0 , t- A) - L(x 0 , t; A) ¥-(x 0 , f; A) (2.22) 

Then using the algebraic relations (2.20) we may formulate on the defect point: 

| In (t); + (L,£o,£; A) Li(x 0 ,t; A) ^(xo, ~L,t; A)), L 2 (x 0 ,L //)} 

= V%(x 0 ,t;\,ii) L 2 (x 0 , t; /i) -L 2 (x 0 ,t;/i) %{x 0 ,t, A, fj) (2.23) 

where ¥ ± are the time components of the Lax pair on the defect point computed from 
the left and right: 

¥^(x 0 ,t; A, fi) = D _1 (xo¥; A) ^l + (L,x 0 ,t; A) r i2 (A - n) L-i (x 0 ,t; A) Th(x 0 , -L,t; A) 
¥ 2 "(x 0 , t; A, n) = D~ 1 (x 0 , t; A) Ti + (L, x 0 , t; A) Li(x 0 , t; A) r 12 {A - n) (x 0 , -L, t; A). 

(2.24) 

Let us recall, as has been discussed in detail in [8], that analyticity conditions imposed 
on the ¥-operators around the defect point, i.e. 

¥±(* 0 ) ^¥±(4) (2.25) 

lead to the necessary sewing conditions on the defect point. These are conditions that 
involve the left and right fields and their derivatives and the defect degrees of freedom. 
Due to these conditions one observes “jumps” of the fields and the derivative across the 
defect point. 

2.2 Integrable systems on the half line 

We shall consider in this section the generalization of the aforementioned algebraic set¬ 
ting for integrable systems on the half line. Before we proceed with the mathematical 
setting in this case it will be instructive to recall the two distinct types of boundary con¬ 
ditions appearing in integrable theories associated to higher rank algebra i.e. gl^. In 
general, depending on the choice of boundary conditions the bulk physical behavior may 
be accordingly affected. In the frame of affine Toda field theories (ATFTs) the boundary 
conditions introduced in [10] are related to the classical twisted Yangian [11], and force a 
soliton to reflect to an anti-soliton, hence the relevant appellation soliton non-preserving 
(SNP) boundary conditions. Naturally another possibility exists, that is boundary con¬ 
ditions that lead to the reflection of a soliton to itself. These boundary conditions are 
associated to classical reflection algebras [16] are known as soliton preserving (SP), and 
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have been extensively investigated in the context of integrable quantum spin chains (see 
e.g. [12] and references therein). Although SP boundary conditions are the obvious ones 
in the framework of integrable lattice models they remained puzzling in the context of 
ATFTs until their thorough analysis in [13]. On the other hand SNP boundary conditions 
were investigated through the Bethe ansatz formulation for the first time in [14], whereas 
higher rank generalizations considered in [15]. 

Let us now describe the algebraic setting associated to integrable systems on the semi¬ 
infinite line. Recall that the modified monodromy matrix for a system on the half line is 
given by [16] 

T(A) = f(A) K( A) T( A) (2.26) 

where we define T in the case one considers the twisted Yangian as: 

f (A) = U T\-\) U (2.27) 

U can be a diagonal or anti-diagonal matrix with particular form depending on the choice 
of the classical r-matrix. We shall later focus on the classical Yangian and consider 
U = antidiag(l, 1,..., 1). In the case one considers the reflection algebra T is defined as 
[16] 

T(A) = T _1 (—A). (2.28) 

We shall focus henceforth on the twisted Yangian case mainly because this case has 
not been really explored in the context of classical continuum field theories such as the 
vector NLS models or the generalized Landau-Lifshitz model. In this case due to (2.26), 
(2.27) it may by easily shown that T satisfies 

{ 71(A), r 2 (/i)} = r 21 (Ai - A 2 ) Ti(A 2 ) T 2 (A 2 ) - Ti(A0 T 2 (A 2 ) r 12 (A x - A 2 ) 

+ 7i(Ai) f i 2 (Ai + A 2 ) 7^(A 2 ) — T 2 (Ai) f i 2 (Ai + A 2 ) 7i(Ai) (2.29) 

where we define in the SNP case (classical twisted Yangian) 

f(A) = Ui r%( A) U\ , (2.30) 

K is a c-number solutions of the twisted Yangian (2.29), 

r 2 i(Ai — A 2 ) Ki(X 2 ) K 2 (\ 2 ) — /^(Ai) K 2 ( A 2 ) r 12 (Ai — A 2 ) 

+Ai(Ai) ri 2 (Ai + A 2 ) A 2 (A 2 ) — A 2 (Ai) fi 2 (Ai + A 2 ) Ai(Ai) = 0 (2.31) 

where here notice that in the sl 2 case the twisted Yangian and reflection algebras coincide. 
Henceforth we shall adopt the twisted Yangian for our computations, extra attention will 
be given to the the gl^, AT > 2 case in this context. In the case associated to the reflection 
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algebra (SP boundary conditions) r '12 = f\ 2 . Although this case is equally interesting it 
involves various technical intricacies and will be left for future investigations. 

In addition to the classical vertex T (2.5) we also introduce 

'b(A) = V f (A). (2.32) 

Moreover, the following relations are required: 

% r 12 (A) V 2 = /(A) ^ V 2 , 

V 2 f 12 (A) Vi = /(A) V 2 V 

Vi V 2 r 12 (A) = /(A) Vi f> 2 (2.33) 

Then in addition to (2.7) the following set of algebraic relations are also satisfied: 

j'bi(Ai), v b 2 (A 2 )| = 'bi(Ai) t 2 i(Ai + A 2 ) V 2 (A 2 ) 
j'bi(Ai), T 2 (A 2 )| = ^/i(Ai) fi 2 (Ai + A 2 ) v b 2 (A 2 ) 

|'Vi(Ai), T 2 (A 2 )} = 'bi(Ai) T 2 (A 2 ) t 2 i(Ai — A 2 ) (2.34) 

where the shifted f-matrix is defined as 

f = f(A)-/(A). (2.35) 


Define the generating function of the integrals of motion as: 

X(A) = V(0, —L, t; A) K{ A) T(0, —L, t; A). (2.36) 

Indeed taking into account the relations (2.7), (2.34) and (2.31) it is shown that 

A), T(A 2 )| =0. (2.37) 


We may now derive the V-operator on the boundary point x = 0 (see also [17] for 
relevant results). The main aim is to formulate the following Poisson commutator taking 
into account the algebraic relations as well as the form of the zero curvature condition at 
x = 0: 


(2.38) 


In (^Ti(0, — L, t ; A) Ai(A) v bi(0, — L, t; A)), U 2 (0, t; 
dV 2 (x,t; A,/i) 


dx 


+ 


V 2 (0,f; A,p), U 2 (0,f;/i) 


where we define the boundary V-operator at x = 0 

V 2 (A,p) = T" 1 (A)(V 1 (A) Ki{ A) r 12 (A-/x) T 1 (A)+T 1 (A) r 12 (A+/x) iV x (A) T X (A)). (2.39) 





Continuity conditions at x = 0 


V(0) ->■ V(0 + ) (2.40) 

lead to suitable boundary conditions on the fields and their derivatives (see detailed 
discussion on [17]). 

Having introduced all the necessary algebraic objects for classical integrable systems 
on the full and half line as well as in the presence of point-like defects we can use a class of 
models associated to the classical (twisted) Yangian r-matrix as a paradigm. Particular 
emphasis will be given on such models with boundary conditions related to higher rank 
algebras. 


3 Paradigm: the classical Yangian & twisted Yangian 

Having at our disposal the general algebraic setting in treating integrable systems via the 
classical version of the vertex algebra we may now consider as a paradigm a whole class 
of integrable models associated to the classical gl^ Yangian with the classical r-matrix 
being [5] 

a r 

r ( A ) = jYl e b ® e H (3- 1 ) 

i,j =1 

in this case /(A) = A 

We shall first consider the model on the full line and demonstrate the simplicity of ob¬ 
taining the integrals of motion. Involution is guaranteed by construction via commutation 
relations (2.10). As already mentioned the fact that the classical vertex operator can be 
immediately obtained from the monodromy matrix provides an elegant way of expressing 
the classical vertex using the standard decomposition of the monodromy matrix [3] 

T(0, -L, t■ A) = (1 + W{ 0, t; A)) e z ^~ L ^ (1 + F(-L, f; A)), 

T*(L, 0, t] A) = (1 + W(L, t; A)) e z(L ’°’ t;A) (l + F( 0, t; A)) (3.2) 

where 

Z = y ^ Zj 6jj , 1 + F = (1 + W) , M = ^ ^ Wij eij, F = y] F^j eij (3.3) 

3 ij 

W and Z may be identified using the fact the monodromy matrix satisfies the first of the 
equations (2.1), and they of course depend on the choice of the model, i.e. the U-operator 
[3]. More precisely, substitution of the decomposition of the type (3.2) for T(x,y, A) in 
the first of the pair of equations (2.1) leads to typical Ricatti equations for W [3]. Solving 
these one can identify W and Z. 
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We assume Schwartz boundary conditions at the end points x = ±L, i.e. the helds 
and their derivatives disappear so that W(±L) = 0. Introduce also the Af dimensional 
column vector u*, with unit on the k th position and zeros elsewhere, also u*. = uj(. Also 
we choose to consider henceforth 


V = uy, V* = uy. (3.4) 

Then the fundamental quantities T, T* with the use of the identities: 

^jk & ill &ij ^0/ &jk U-j, h/,. Cjj A A'hi. U , Uy Sjj (3.5) 

become 

T(0, —L, t; A) = e z ^~ L ^ (i W + ^ W jN { 0, U A) Uj ) 

T*(L,0,f;A) = e z ^ L ’°^(u A , + ^F A , J (0,t;A) fq). (3.6) 


Note that in general we could have used the standard inverse methodology to identify the 
the rations |^, where obviously and equivalently to (3.6) 


AT 


AT 




U) 


Ui 


(3.7) 


3 = 1 


3=1 


Indeed, recalling that T satishes the first of the equations (2.1) and expressing the U 
operator as 

U(x, t; A) = Ujj(x, £; A) , (3.8) 




one obtains 

d_ 

dx - 


e Z M {x, L,t) u ^ +e Zjs(x, L,t)^2w jAr (x,t)Uj = Y, Uiujx, t)uy+ Y, U ij(x, t)W jA f(x, t)Uj 

j¥=Af i j^Af,i 

(3.9) 

which leads to the typical Ricatti type equations for already mentioned above 
dWkJ ^ x ' t) + (lUyy^f) + ^D^(i,t)W>(i,t)]%(i) = 


3+ 1 


u kAr(x, t) + ^2 t)Wjtf(x, t ) 


3+A[ 


= U»( I( () + E U w (i, f)»>(x, i), 


(3.10) 


3+AI 


solution of the equations above leads to the identification of Zy, W^y. 
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The advantage here is the use of the decomposition of the monodromy matrix (3.2), 
which provides the simple formula below for the generating function of the integrals of 
motion: 

ln(T(L, -L, t; A)) = Z^(L, —L, t; A) (3.11) 

with apparent analyticity conditions at x — 0, this will become more transparent when 
discussing such models in the presence of point-like defects later in the text. It will be 
also instructive to derive the V-operator (2.17). Taking into account (2.17), (3.1), (3.6) 
we conclude that 

V(x,f; A,/i) = - ^ + FujW e A/j + ^^ Wjj\f( A) e^-\- Fjj *(A) Wjx{\) 

^ 3 i-3'3 i 

(3.12) 

which coincides with the generic expression for V in the context of classical Yangian [3, 8]. 
Substitution of the W % j and Z^j quantities on the expressions for the generating function 
and the V-operator provides explicit expressions at each order for a particular model 
under study (see e.g. [3, 8] for explicit expressions). 

3.1 Defects 

The main aim in this subsection is to identify explicit expressions of the generating func¬ 
tion and hence the local integrals of motion as well the time component of the Lax pair 
using the classical vertex operators in the presence of a local defect for the classical Yan¬ 
gian case. 

Recall the expression for the generating function (2.18), it is then straightforward to 
rewrite it via (2.19), (3.6), as 

ln(D(A)) = Z^(L, x 0 , t; A) + Z^{x 0 , —L, t; A) + lnX(xo, t ; A) (3.13) 

the plus and minus superscripts in the expression above refer to the left right bulk theories, 
and define 

X(*o, t; A) = LawW+^T Ffo (A) L^(A)+Y L W (A) W^(A)+£ F&( A) L«(A) Wr„( A) 

3 jjM 3+N 

(3.14) 

Similarly the time components of the Lax pair on the defect point (2.24), take the 
explicit expressions after substituting (3.1), (2.19), (3.6) in (2.24) 

V (A, /i) — — - ^ 'y T y T y F^-^Luj^eij + Y F^W^hue,,) 

^ 3 Vlj *7■WjV 
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^ (^j A 4 ) ^ _ f ^ 1 ^Mj e Mj ~b ^ ^ A ^ 1 -^W^ja e .A^ ~b /* 1 

^ j 'Y-Yi i,j i^AT,j,k 

The expressions of the generating function and the defect V-operators coincide as expected 
with the ones of earlier works on models with underlying classical Yangian [8]. Explicit 
expressions associated e.g. to the NLS model and a models in the presence of defects can 
be found in [8], and they are not repeated here for brevity. The interested reader however 
is referred to [8] for more details. 

3.2 The system on the half line 

We shall now derive explicit expressions for any system associated to the classical twisted 
Yangian such as the vector NLS model or the isotropic Landau-Lifshitz model or principal 
chiral models. Note that a systematic study at the level of classical field theories is still 
missing although there are some preliminary results in [18], but mainly focused on the 
discrete analogues of these models. This is in fact the first time that this problem is 
directly addressed at the level of continuum classical systems. 

Let us recall the generating function of the integrals of motion in this case given in 
(2.36), also in addition to the explicit expression for T (3.6) one needs the expression for 
T. Indeed choosing to consider V = uy we obtain 

T(0, -L, t; A) = ^ + ^ ^(0, f; A)u,) (3.15) 

m 

where: Zy(A) = Z[y(—A), W i3 { A) = W l3 (—X), and we also define the conjugate index j = 

M — j T 1. Also, the c-number /i-matrix can be generically expressed as: K = JT . 
then the generating function becomes 

‘1(A) = e z *W+z*W y(A) 

Y(A ) = K 1Ar +Y,Ki j W j „+Y,K j xW- j „+ J2 W ]M K p W lM . (3.16) 

i+Af j^AT jYP'YA/' 

Due to the generic expression for the boundary V-operator (3.6), (3.15) and 
Y(0,f;A,/r) = 7 [£(*■+£ KwWij^je^fi + (-^i iWiM + ^2 KjiWjAfWiN^eu j 

^ 1 l 1 i+M,l V 1 

— — [ ^2 (-^Tv + ^2 e M + ^2 (K-lM'WjAf + ^2 K-IiWjj\rWij2) e ]i 

+ ^ i j&r m 

(3.17) 
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It is worth presenting here, in addition to the generic expressions (3.16), (3.17), specific 
conserved quantities associated to the vector NLS model. In [18] results on the boundary 
vector NLS model, but using its discrete counterpart and then taking a suitable continuum 
limit. Let us now focus on the vector NLS model; the corresponding U-operator 

U(A) = U 0 + AUi (3.18) 

where U; are AT x AT matrices 

^ AT-1 

Ui = 2^( e ' J1 ~ e AOV^ 

3 =1 
Af-1 

U 0 = y/R ^2 (i’jejM + V’je.vj) (3-19) 

3 = 1 


Then from the Riccati equations one obtains (see also [18]) 

W^f = . . . (3.20) 

and 

f B 

Ztf(B,A,t;A) = iL8 n - 1 + y/Kj dx ipj(x, t) t; A) (3.21) 

Given the expression for the local charges ln(T) (3.16) we expand in powers if j- and 
the first non trivial conserved quantity is the modified momentum 

AT-l o _ _ AT-i _ M-i 

J (2 ) dx yPj( x ) Vj( x ) “ ^j( x ) ^A x )) + K X] ^(°) + « ^2 ^i(°) ^(°)- 

j =1 j =1 3 =1 

(3.22) 

This result confirms earlier findings on the vector NLS model, which were obtained from 
the study of the discrete version of the NLS model [18]. Application of this setting in 
other physical systems associate to higher rank algebra such as affine Toda field theories 
or higher rank a models would certainly provide significant physical results, but these 
issues will be addressed in forthcoming investigations. 


4 Discussion 

We have formulated a classical variation of the quantum vertex algebra, and we have 
shown that the auxiliary linear function is basically a classical vertex operator. Based 
on these notions we considered the problem of identifying the local integrals of motion 
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and time components of the corresponding Lax pairs. Particular emphasis was given in 
systems in presence of defects and non-trivial integrable boundary conditions. It turns out 
that in these cases and in particular in the boundary case this formulation is technically 
considerably simpler and provides results directly at the continuum level even in situations 
that this has not been possible in a straightforward manner such as e.g. in the study of 
the vector NLS with twisted Yangian symmetry. 

One of course has to note that the choice of the vectors V, V* implemented at the 
ends of the system modifies or even perhaps restricts the behavior of the system under 
study. It is worth pointing out however that complications arising due to the asymptotic 
behavior of the monodromy matrix in systems associated to high rank algebras on the 
half line can be efficiently dealt with in this framework, basically due to the restrictions 
imposed from the vectors V, V*. In fact, it would be illuminating and would provide 
further physical insight to implement this frame in prototype classical systems such as 
affine Toda field theories with integrable boundary conditions, however this issue will be 
addressed in a separate publication. 

At the quantum level and in the context of quantum inverse scattering method in par¬ 
ticular one builds the one dimensional lattice system using the “bare” A-matrix. Then 
the solution of the Bethe ansatz equations (BAE) in the thermodynamic limit serves as 
a renormalization process, so one may for instance derive the physical scattering infor¬ 
mation solving the BAE in the thermodynamic limit. Formulation of a bare version of 
vertex operators and the corresponding vertex algebra using the bare A-matrix could be 
the first step towards reconciling the Bethe anaszt formalism with the vertex operator 
construction by the Japanese group [21], The light cone lattice approach by Destri and 
de Vega [22] could serve as a possible bridge between these two methodologies. This is a 
very interesting direction to pursue and will be left for future investigations. 
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